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BIFURCATIONS OF CRITICAL ORBITS OF INVARIANT 
POTENTIALS WITH APPLICATIONS TO BIFURCATIONS OF 
CENTRAL CONFIGURATIONS OF THE N-BODY PROBLEM 

MARTA KOWALCZYK 


Abstract. In this article we study topological bifurcations of critical orbits of equivari- 
ant gradient equations. We give necessary and sufficient conditions for the existence of 
global bifurcations of solutions of these equations. Moreover, we apply these abstract 
results to the study of bifurcations of new families of planar central configurations of the 
N-body problem. 


1. Introduction 

The most important questions of celestial mechanics are the ones concerning central con¬ 
figurations. They have a long history dating back to the 18th century and have been 
studied by various mathematicians, including Euler and Lagrange. The first known cen¬ 
tral configurations were the three collinear classes discovered in 1767 by Euler (see [11]) 
and the two classes with masses located at the vertices of an equilateral triangle found in 
1772 by Lagrange (see [18]). Later it was proved by Wintner that the 3 collinear and 2 tri¬ 
angular configurations are the only possible central configurations in the 3-body problem 
(see [38]). In general, the question about the finiteness of number of central configurations 
has been included in the famous list of eighteen problems of the 21st century by Smale 
([35]) and it is the following: 

“Given positive real numbers m i,..., m n as the masses in the n-body prob¬ 
lem of celestial mechanics, is the number of relative equilibria finite?” 

As a matter of fact, this open question was already formulated by Wintner in his book 
([38]) in 1941. In the 4- and 5-body problem the answer to Smale’s question ([35]) has been 
given by Hampton and Moeckel ([16]), and by Kaloshin and Albouy ([1]), respectively. 
Although as for now the question remains unanswered for 6 and more bodies, there were 
given some estimations for the number of central configurations. For instance, Merkel in 
[22] and Pacella in [27], using Morse theory, determined a lower bound for the number of 
spatial non-collinear central configurations (for a lower bound of the number of planar non- 
collinear central configurations see [28]). For a deeper discussion of central configurations 
we refer the reader to [26]. 
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Classifying central configurations is a very difficult problem, therefore an answer is 
sought even with some simplifications, that is imposing restrictions on the masses or 
considering only highly symmetrical solutions. One of the many people to go in this 
direction was Palmore, who considered for n — 4, 5 the configuration of n — 1 bodies of 
mass 1 at the vertices of a regular (n — l)-gon and n-th body of arbitrary mass m at 
its centroid ([29]). He showed that for n — 4 and for n — 5 there exist unique values of 
the mass parameter m for which this central configuration becomes degenerate. Meyer 
and Schmidt in 1988 reproduced these results and additionally showed that a bifurcation 
occurs ([24]). Since then a lot more central configurations have been discovered in spirit 
of these mentioned restrictions. Examples of central configurations which are considered 
in our paper are the following: 

• a planar configuration of two nested squares found in 2013 by A. C. Fernandes, L. 

F. Mello and M. M. da Silva in [12], 

• a planar rosette configuration of 13 bodies studied in 2006 by J. Lei and M. 

Santoprete in [19], and in 2004 by M. Sekiguchi in [34], 

There are numerous reasons for studying central configurations. There exists an im¬ 
portant class of solutions of the n-body problem - the so-called homographic solutions 
([38]) - which is directly linked to the concept of central configurations. A solution is 
called homographic if at every time the configuration of the bodies remains similar. As a 
matter of fact, this can only happen for central configurations. The Euler and Lagrange 
solutions, which we have already mentioned, are the most famous examples. Thus cen¬ 
tral configurations generate the unique explicit solutions of the iV-body problem known 
until now. In particular, planar central configurations also give rise to families of periodic 
solutions. As a special case, the homographic solutions for which the configuration just 
rotates rigidly at constant angular speed are known as relative equilibria. Moreover, cen¬ 
tral configurations are important for instance in analysing total collision orbits, namely 
a certain change of coordinates in the neighbourhood of a total collision point turns the 
colliding particles into a central configuration in the limit as t approaches the collision 
time. This phenomenon was shown initially for the complete collapse of all bodies at 
the centre of mass by Wintner in [38] and was later generalised by Saari in the study of 
the parabolic escapes of the particles. The same holds for expanding subsystems, which 
also asymptotically tend to a central configuration (readers interested in these two phe¬ 
nomena should see [33] for details and more references). These examples demonstrate 
that knowledge of central configurations gives important insight to the dynamics near a 
total collision and asymptotic behaviour of the universe. Another aspect motivating the 
search for these special solutions is the fact that the hypersurfaces of constant energy and 
angular momentum change their topology exactly at the energy level sets which contain 
central configurations ([36]). 

Central configurations are invariant under homotheties and rotations - the first claim 
comes from the homogeneity of the potential while the second one stems from the fact 
that the potential depends only on the mutual distances between the particles and not 
on its positions. Additionally, in the planar case the action of the symmetry group, in 
this case SO( 2), is free, so the quotient space Vl/SO( 2) of the action of the group SO( 2) 
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on the configuration space is a manifold. Because of this a natural thing to consider is 
a new potential defined on the quotient space Q/SO( 2), which has been widely used in 
the literature ([19, 23, 24, 25]). In this context central configurations have been viewed 
not as critical orbits, but as critical points of a “quotient” potential. On the other hand, 
in the spatial case the action of the symmetry group 50(3) is not free, therefore the 
quotient space 0/50(3) of the action of the group 50(3) on the configuration space O is 
not a manifold (there are two isotropy groups appearing in the configuration space), so 
the use of ordinary theory is not justifiable (without any differential structure). Thus in 
this case we cannot use the approach from the planar one. For example, the equivariant 
Morse theory has been presented to the study of spatial central configurations in [27]. In 
general, a different approach which does not employ the concept of quotient space has 
also been successfully used ([14, 20, 30]). Apart from the symmetries, Newton’s equations 
are of gradient nature, therefore it seems reasonable to apply methods which benefit from 
this. 

In this paper we use the second approach to study classes of planar central configura¬ 
tions, which are treated as SO (2)-orbits of critical points of 50 (2)-invariant potential. 
Among the machinery which is designed to work in the equivariant context making use 
both of the symmetry and the gradient structure appearing in the problem and which 
we employ in our paper are equivariant versions of classical topological invariants such as 
the degree for equivariant gradient maps (see [15] and [32]) and the equivariant Conley 
index (see [4], [13] and [15]). This approach using equivariant theory has for instance been 
used in [20] and [30]. For general treatment of the degree for equivariant maps we refer 
the reader to [2] and [3]. In this article we apply those topological tools in equivariant 
bifurcation theory to produce new classes of planar central configurations from known 
families by obtaining global bifurcations of families of planar central configurations under 
a change of the degree for equivariant gradient maps. The same result has been reached 
for local bifurcations under a change of the equivariant Conley Index ([37]). In the planar 
case because we have exactly one orbit type on the configuration space, the results ob¬ 
tained by using the equivariant Conley Index or the degree for equivariant gradient maps 
are equivalent to those which one could get considering the quotient space and applying 
ordinary Conley Index or the Brouwer degree. However, in the spatial case the approach 
via the quotient space is not possible. In general, the obtained theoretical results are 
stronger than these we need to study of the planar ./V-body problem. 

The aim of this paper is twofold. On the one hand, we give new results about equivariant 
bifurcation theory (see Theorems 3.2, 3.3 and 3.4). On the other hand, we apply those 
theoretical results to the study of planar central configurations. In the next article we 
shall present an application of our abstract results to the spatial ./V-body problem and 
give necessary and sufficient conditions for the existence of local and global bifurcations of 
spatial central configurations. Thus the main goal of this article is to create a homogeneous 
theory which can be used for both planar and spatial At-body problems. In the planar 
case, as we mentioned before, we can consider the orbit space of the action of the group 
50(2) on the configuration space and then use both ordinary Morse theory or the Brouwer 
degree and the equivariant approach. However, in the spatial case we cannot work on the 
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quotient space because the action of the group 50(3) on the configuration space is not 
free. 

Notice that we can consider the possibility of working on subsets of the full config¬ 
uration space which are invariant for the gradient flow and seek bifurcations in these 
subsets, but it is possible that we can get no bifurcations while there are ones in the full 
configuration space. This phenomenon happened for instance in [24], that is bifurcations 
of less symmetrical families from highly symmetrical families of central configurations 
occur. Meyer and Schmidt treated the highly symmetrical family of equilateral triangle 
with fourth body at the centroid and proved that another families which are less symmet¬ 
rical (that is the families of isosceles triangles with fourth body near the centroid and on 
the line of symmetry of the triangle) bifurcate from this family. Similarly, these authors 
taking the family of square with fifth body at the centroid proved bifurcations of families 
of kites and isosceles trapezoids (for more examples of this phenomenon see [25]). In our 
case of the family of rosette configuration of 13 bodies, likewise, we have the existence of 
bifurcations of central configurations (not necessarily configurations of rosette type) from 
the rosette family (see Theorem 4.2) and it transpires that central configurations which 
bifurcate are not of rosette type. However, we do not know the shape of the bifurcating 
families. One can compute that there are no values of the mass parameter for which the 
rosette family is degenerate. Therefore the obtained results are stronger than these which 
could be get by taking into account the invariant subsets of the configuration space. In 
the case of the family of two nested squares we also get stronger results not considering 
invariant subsets. 

We stress the fact that the algebraic structure of the potential is not important from the 
point of view of central configurations, the only things that matter are the rotation and 
scaling invariance. This makes the methods that we use applicable to mathematical mod¬ 
els of physical problems involving interactions between multiple bodies, whose behaviour 
is not necessarily governed by Newton’s gravity laws. Such situation arises for instance in 
molecular dynamics, where intermolecular relations are modelled using a broad spectrum 
of potentials. The famous Lennard-Jones potential, which approximately describes the 
behaviour of two neutral atoms or molecules, was studied in this context in [6] and [7]. 

This paper is organised as follows: Section 2 contains an introduction to the equivariant 
setting along with the definition and properties of a topological degree suitable in this 
context and lemmas used in the proofs of our theorems. In Section 3 we investigate the 
G-orbits of solutions of the equation 

V v p(v,p) = 0 (1.1) 

under the action of G, strictly speaking we formulate abstract theorems giving necessary 
and sufficient conditions for the existence of bifurcations in the vicinity of a given family 
of solutions. In Theorem 3.1 we state that only a degenerate critical G-orbit of solutions 
of equation (1.1) can be a G-orbit of local bifurcation. The sufficient condition is given 
by Theorem 3.2, its assumptions, which include a change of the degree for G-equivariant 
gradient maps computed at critical G-orbits of p, imply the existence of a global bifurca¬ 
tion. We also provide simple conditions in Theorem 3.3 and Remark 3.1, which allow us 
to verify inequality of degrees. Moreover, we prove a theorem in the spirit of Rabinowitz’s 
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famous alternative, see Theorem 3.4. In Section 4 we apply our abstract results to the 
planar iV-body problem, namely we consider known families of central configurations - 
the two squares family (4.4) and the rosette family (4.5) - and prove the existence of bi¬ 
furcations of new classes of central configurations from these families. We stress that our 
method is applicable not only in this specific example, but is rather a general framework, 
which is additionally easily implementable on any computer algebra system, which allows 
symbolic computations. Notice that in this paper we consider topological bifurcations 
(local and global ones, see Definitions 3.1 and 3.2) rather than bifurcations in the sense 
of a change in the number of solutions. In other words, in the case of local one, we 
consider known family of solutions of equation (1.1) (so-called trivial solutions) and seek 
non-trivial ones nearby. In the case of global one, we seek connected sets of non-trivial 
solutions nearby the trivial ones which additionally satisfy Rabinowitz type alternative 
(see condition (3.3)). Those two definitions, the bifurcation in the sense of a change in 
the number of solutions and the topological bifurcation, are independent, that is the first 
one can occur while the second one does not occur and inversely. Throughout this paper 
we will write about topological bifurcations briefly bifurcations. 

2. Preliminaries 

In this section we review some classical facts on equivariant topology. The material is 
fairly standard and well known (for details see for example [10] and [17]). Throughout 
this paper we will remain in the real setting, that is all vector spaces will be considered 
over M. and all matrices (treated as elements of Lie groups) will have real entries. Also, 
from now on G stands for a compact Lie group and V = (M n ,?) is a finite-dimensional, 
real, orthogonal (S'-representation, that is a pair consisting of the Euclidean n-dimensional 
vector space and a continuous homomorphism g : G —» O(n), where 0(n ) denotes the 
group of orthogonal matrices. By v G V we mean v G M n . In addition to that by SO(n ) 
we will understand the group of special orthogonal matrices. The linear action of G on 
M n is given by £ : G x M n —y M n , where £(g, v) = g(g)v. We write gv for short. Moreover, 
M denotes the one-dimensional, trivial representation of the group G, that is M = (M, 1), 
where 1(g) = 1 for any g G G. Then the linear action of G on M n x 1 is given by the 
Cartesian product of <, and 1 

G x (R n x R) 3 (g, (v,p)) ^ (gv,p) G M n x R. 

A subset Q C V is called (S'-invariant if for any g G G and v E Q we have gv G D. 
Let sub(G ) be the set of closed subgroups of G. Two subgroups H, K G sub(G ) are called 
conjugate if there exists g G G such that H = g~ l Kg. Conjugacy is an equivalence relation, 
the conjugacy class of H is denoted by (H) and it is called the orbit type of H. Let sub[G] 
be the set of conjugacy classes of closed subgroups of G. Similarly, H is subconjugate to 
K if H is conjugate to a subgroup of K, denote by (H) < (. K ). Subconjugacy defines a 
partial order in sub[G\. Additionally, we write (H) < (.K ) if ( H) < (K) and (H) ^ (. K ). 

If v G V then G v = {g G G : gv = u} G sub(G) is the isotropy group of v. For each 
v G V the set G(v) = {gv : g G G} C V is called the G'-orbit through v. The isotropy 
groups of points on the same G-orbit are conjugate subgroups of G. We will also denote 
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by V G the set of fixed points of the G-action on V, that is V G = {v G V : G v — G } = 
{v G V : gv = v \/g G G}. Moreover, for given H G sub(G) and an open, G- invariant 
subset O C ¥ we have Ul<(H) = {» G O : (G v ) < (//)} and £1(h) = {v G Q : (G v ) = (//)}. 
Let H G sub(G ) and W be an orthogonal //-representation. Now define an action of H 
on the product G x W by the formula (h, (g,w)) (->• ( gh _1 , hw ) and let G Xh W denote 
the space of H -orbits of this action. 

Fix k,l G NU { 00 } and an open, G-invariant subset hi C V. A map p : hi x M —» M 
of class C k is said to be a G-invariant C k -map if < p(gv,p ) = < p{y,p ) for any g E G and 
(v,p) G hi x M. The set of G-invariant G fc -maps is denoted by Cq(G, x M,M). A map 
if} : hi x M —y V of class C l is said to be a G-equivariant C l -map if ip(gv,p) = g^(v,p) 
for any g G G and (v, p) G 0 x R. The set of G-equivariant G^-maps is denoted by 
C 1 q{VL x M, V). For any G-invariant G fc -map p the gradient of p with respect to the first 
coordinate denoted by V v p is G-equivariant G fc-1 -map. 

Let OcVbe open and G-invariant and £x vq G hi, where by H we denote G Vo . Notice 
that Q(h) is a G-invariant submanifold of hi and G(v 0 ) is a G-invariant submanifold of 
f 1(H), therefore we obtain the following orthogonal direct sum decomposition 

T Vo V = T Vo Q = (T Vq Q(h)) © (TvqQ^h )) 1 = (' T Vo G(vo )) © (T Vo Q(h) QT Vq G(v 0 )) © ( T Vo fl^ H )) ± , 

( 2 . 1 ) 

where T vo V is the tangent space to V at Vq. On the other hand, we have 

T vo V = T vo n = (T V0 G(v 0 )) ©W ff ® (W^) 1 , (2.2) 

where W = (T„ 0 G(uo))" L is an orthogonal //-representation. Because (Gx^B; o (W))(h) = 
G/H x B £ vo (W H ), we get that W H C T Vq Q( H ), where B £ vo { W) denotes the e-ball centred 
at V(j in W. Now assume that p G Gg(h2,M) and Vq G (V(/9)^ 1 (0). The following lemma 
gives the decomposition of the Hessian V 2 p of the potential p with respect to the sum 
decompositions given by the formulas (2.1) or (2.2) and has been proved in [15]. 


Lemma 2.1. Under the above assumptions, the Hessian 


V 2 p(v 0 ) : T V0 Y = 


Ty 0 G{v 0 ) 

Ty 0 G(v 0 ) 

© 

© 

W H -> Ty 0 V = 

W H 

© 

© 

(W H ) X 

(W H ) X 


is of the form 


"00 0 
VV(u 0 )= 0 B(v 0 ) 0 

0 0 C(v 0 ) 


(2.3) 


According to Lemma 2.1, we have dimker V 2 p(vo) > dimG(uo). We will thus call 
a critical G-orbit G(v 0 ) degenerate if the strict inequality holds and non-degenerate if 
dimker V 2 </j(uo) = dimG(u 0 ) (which is equivalent to B(v 0 ) and C(v 0 ) being non-singular). 
Note that non-degenerate critical G-orbits are isolated, that is there exists a G-invariant 
neighbourhood 0 C V of G-orbit G(v 0 ) such that (V(/9) _1 (0) D 0 = G(v 0 ). Moreover, a 
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non-degenerate critical G-orbit G(v o) is called special if m~(C(v o)) = 0, where by the 
symbol m~(C(v o)) we denote the Morse index, that is the number of negative eigenvalues 
(counting multiplicities) of the matrix C(vq). 

Regarding planar central configurations, there is exactly one orbit type (H) in the 
configuration space hi, that is hi = Q(h) ((H) = ({e})), so to illustrate this situation we 
give the following remark and example: 

Remark 2.1. Under the assumptions of Lemma 2.1 with hi = h 1(h), the Hessian 

T V0 G(v o) T V0 G(v o) 

VV(n 0 ) :©—>■© 

T Vq Q(h) © T Vo G(v o) T Vq Q(h) © T Vq G(vq) 

is of the form 

V 2 (p(v o) = 

where B(v 0 ) is from the formula (2.3). 

Example 2.1. Let G = SO(2) and V = (M 2 ,^), where q : SO( 2) —» 0(2) is given by 
q(g) = g. Assume that hi = V\{0}, then for any v G hi we have SO(2) v = {e}, so 
hi = hl(( e |). Fix ip G 0| 0(2) (hl,M) given by (p(v) = ip(\v\ 2 ), where ip G G°°(M,M) is dehned 
by ip(t) = t(t — 1) and by the symbol |n| we denote the standard norm of an element v. 
Then V<p(v o) = 0 if and only if |n 0 | 2 = p, so we get that for v 0 = (0, ^) the Hessian 

T V0 SO( 2)(v 0 ) T V0 SO( 2){v 0 ) 

X 2 <p(v o) : © —y ffi 

^ 0 ^({e}) © T Vo SO(2)(v 0 ) T„ 0 hl(( e p © T Vo SO(2)(v 0 ) 

is of the form 

V 2 <p(v 0 ) = 

Note that if there is only one orbit type in hi then the quotient space hl/G of the action 
of the group G on the open and G-invariant set hi is a manifold, so for any G-invariant 
G 2 -map we can consider a "quotient" potential dehned on hl/G. In the following lemma 
we present the decomposition of the Hessian of the "quotient" potential and it transpires 
that the ordinary Morse index of the "quotient" potential is equal to the Morse index of 
matrix B(v o) (see the formula (2.3)). 

Lemma 2.2. Let hi C V be open and G-invariant. Assume that hi = LI(h) for some (H) G 
sub[G] and fix p G Gg(hl,M). Then the potential ip : hl/G —>■ M given by ip(G(vfi)) = <p(vq) 
is a map of class G 2 and V 2 ip(G(v o)) : T Vo Ll/G —> T„ 0 hl/G is of the form . V 2 ip(G(v o)) = 
[-B(no)], where B(v o) is from the formula (2.3). 

Let BfiG) denote the class of finite, pointed G-G W -complexes (for the definition of 
G-GhU-complex see [10]). The G-homotopy type of X G J-/(G) is denoted by [A"]. Let 
F be a free abelian group generated by G-homotopy types of finite, pointed G-CW- 
complexes and let be a subgroup of F generated by elements [A] — [A"] + [A /A \, where 
A, X G A*(G) and A C X. Put U(G) = F/N and let xg(X) be the class of an element 


0 0 
0 4 


0 0 
0 B ( v 0 ) 


(2.4) 
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[A"] G F in 27(G). If X is a G-CW- complex without base point we put xg(X) = xg(X + ), 
where X + = X U {a}. The ring U(G) is called the Euler ring of G ( for the definition of 
27(G) see [9] and [10]). The following theorem can be found in [10]. 

Theorem 2.1. The group ([27(G), +) is the free abelian group with basis xg(G/H + ) for 
(77) G sub[G\. Moreover, if X G 77(G) and [Jf, =0 {(A;, (77^)) : j = 1, • • • ,q(k)} is a type 
of the cell decomposition of X, then 

Xg{X) = X (£(-l )G(k,(H))Yxa(G/H+), 

( H)£sub[G] \k=0 J 

where u(k , (77)) is the number of cells of dimension k and of orbit type ( H ). 

Let 12 C V be open and (S'-invariant. A function p G G g (V, M) is said to be 12- 
admissible if (V</ 3 ) - 1 ( 0 ) D <912 = 0. An 12-admissible map p G G g (V, M) is called a special 
Morse function if for any vq G (V(/ 3 ) _ 1 ( 0 ) fl 12 a critical G-orbit G(v o) is special. 

Let 12 C V be open, G-invariant and bounded. For an 12-admissible function p G 
Cq(Y, M) there is a notion of degree for equivariant gradient maps V G -deg(V</?, 12) G 
U(G). This invariant has been introduced by G^ba in [15] with coordinates as follows 

V G -deg(Vv?, 12) = ^ V G -deg w (Vy 9 ,12) • Xg (G/77 + ) G U(G). 

(H)e^ub[G] 

To make this article more readable we give a sketch of the definition of the degree. First, it 
is proved that there exists a special Morse function if G Cq(Y, M) which is G-homotopic 
to the map p, that is there exists H G G g (V x [0,1],M) being 12-admissible for each 
t G [0,1] such that H(-, 0) = if and H(-,l) = p. Additionally, because non-degenerate 
critical G-orbits are isolated and 12 is bounded, we get that (VV ; ) _ 1 (0) fl 12 consists of a 
finite number of distinct critical G-orbits, that is (V - ?/ ; ) _ 1 (0) fl 12 = G(tq) U ... U G(vi), 
where G{vf) fl G{vf) = 0 for i j. Let us remind that for any special critical G-orbit 
G{vf) we have m~{Cfvf)) = 0, so next we allocate (—l) m ^Gd) to each critical G-orbit 
G[vi ) and define the degree for equivariant gradient maps V G -deg(V(/?, 12) G U(G) by 

V G -deg (Vp, 12) = V G -deg(V^, 12) E (-l) m " (B(,;i)) I -Xg(G/H + ) g U(G). 

(H)esub[G] \(G Vi )=(H) ) 

The definition of V G -deg(V(p, 12) does not depend on the choice of the special Morse 
function. For details of the definition and properties of the degree for equivariant gradient 
maps see [15] and [32], 

In the following theorem we formulate some basic properties of this degree. 

Theorem 2.2. Fix an 12- admissible function p G G g (V, M), where 12 is an open, bounded 
and G-invariant set. Then the degree V G -deg(V(/ 2 ,12) G U (G) has the following properties: 

(1) [Existence] if V G -deg(Vy>, 12) / 0 G U{G), then (V(/ 9 ) _ 1 ( 0 ) fl 12 7 ^ 0, 

(2) [Additivity] if 12 G 122 are open, bounded, G-invariant and disjoint sets such that 
12 = 12i U 122 , then 

V G -deg (Vp, 12) = V G -deg(V<^, 12 x ) + V G -deg(V<^, 12 2 ), 
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(3) [Homotopy invariance] if $ G Cq(Y x [0, 1], M) is Q-admissible for each t G [0, 1], 
then 

V G -deg(V„$(-,0),fi) = V G -deg(V„$(-,l),fi), 

(4) [Excision] if ft i C 0 is an open, G-invariant set such that (Vy?) _ 1 (0) D fl C fli, 
then V G -deg(V</ 9 , Q) = V G -deg(Vy?, fli). 


In the non-equivariant case there are some invariants for continuous maps, for instance 
the Brouwer degree. This invariant is sometimes too weak to distinguish homotopy classes 
of two equivariant gradient maps. On the other hand, they can be distinguished by the 
degree for equivariant gradient maps. In the example below we define two equivariant 
gradient maps whose the Brouwer degrees are the same while they are distinguished by 
the degree for 50 (2)-equi variant gradient maps. 


Example 2.2. Let G = 50(2) and V = (M 2 ,^), where q : 50(2) —» 0(2) is given by 
q(g) = g. Notice that for any v G V we have 



{e}, v ± 0 
50(2), u = 0 


and £>q(V) C V is an 50(2)-invariant set. Consider two maps ip± G Cg 0 ^(B[(V),R) 
given by the formulas (p±(v) = ±||u| 2 . According to Lemma 2.1, for 0 G (V 93 ±) _ 1 ( 0 ) we 
have 


vV( 0 ) = [ c + ( 0 ) ] 


and VV-(O) = [ C_( 0 ) ' 


-1 0 
0 -1 


Then the Brouwer degrees of the maps V</?± at 0 in the set B[(V) are equal and 


deg B (V 93 +, (V), 0) = deg B (V</?_, .Bq (V), 0) = +1 G Z. 


Because tp + is a special Morse function and it has only one critical 50(2)-orbit 5O(2)(0), 
we get 


V c-deg( V (/?+, Bp (V)) = (-l)’“-< fl +<»)) Xso , 2) (SO(2)/SO(2)+) = 

= (-1)°xso ( 2)(SO(2)/SO(2)+) = 1 e [2(50(2)). 

Next, let H G Cg 0 ^(B^(V) x [0,1], M) be given by H(v,t ) = (1 — t'y(\v\))(p_(v) + 
f 7 (|u|)</? + (u), where 7 G 0°°(M,M) is such that 7 |(_ 00 i i] = 1, 7|[i i + 00 ) = 0 and 7|(i,i) is a 

smooth, decreasing map. Notice that H(-, 0) = </?_, H(-, 1) = (1—7( | • I)) V 7 — (')+7( I * I) V p -f(‘) 
and H is B ( \ (V)-admissible for each t G [0,1], so the homotopy invariance property of the 
degree (Theorem 2.2.(3)) implies that V G -deg(V</?_, B[(Y)) = V G -deg1), B[(Y)). 
One can compute that the map H(-, 1) is a special Morse function and it has exactly two 
critical 50(2)-orbit 5O(2)(0) and 50(2)(u 0 ), where v 0 = (0 ,y 0 ) and \ < y 0 < 1. Accord¬ 
ing to Lemma 2.1, for 0 G (V„iL(-, 1)) _1 (0) and v 0 G (V v H(-, 1)) _1 (0) we have 

0 0 
0 B(v 0 ) 


V 2 5T(0, 1 ) = VV(0) and V 2 5T(u 0 ,1) 
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Then 

V G -deg(V„^(-,l), J B 0 1 (V)) = 

= (- l) m ^ B ^ Xso{2) {SO{2)/{e} + ) + (—l) m ^ B+ ^xso{2)(SO(2)/ 50(2) + ) = 

= ±xso{2){SO(2) / {e} + ) + X so { 2){SO(2)/SO(2 ) + ) G U{SO{2)). 

Now we consider equivariant gradient maps whose the Brouwer degree vanishes and the 
degree for equivariant gradient maps can be nontrivial, see the following example: 

Example 2.3. Let G = 50(2) and O C V be an open, SO (2)-invariant and bounded 
subset of an orthogonal SO (2)-represent at ion V such that Q SO( - 2 ' 1 = 0. Then for any 
Q-admissible map p G 0^.(fl,M) we have degB(V</?, fl, 0) = 0 G Z. 

Below we formulate lemmas which are standard in equivariant topology and will be 
useful for us in further investigations. 

Lemma 2.3. Fix an element g G G. Then the map : V —> V given by the formula 
7 g (v) = gv is a homeomorphism. In particular, the set GU = U g ^GTg{U) C Y is open 
and G-invariant for any open subset U C V. 

Lemma 2.4. Let p G Cq(V, M) and assume that U C V. Then (Vv?) _1 (0) (1 GU = 0 if 

(v<^)- 1 (o)nf/ = 0. 

The following theorem can be found in [5]. 

Theorem 2.3. Let K be a compact space and A, B C K be closed, disjoint sets such 
that there is no connected set S C K such that 5 HAL ^ 0 and S DB ^ 0. Then there exist 
compact sets Ka, Kb C K such that A C Ka, B C Kb, KaOKb = 0 and KaOKb = K. 

Lemma 2.5. Fix a map p G Cq(Y, R) and v 0 G V such that Vtp(vo) = 0. Then for any 
v G GBj, o (yV) we have (G v ) < (G Vo ), where W = (T„ 0 G(uo))" L - 

The following lemma can be found in [21]. 

Lemma 2.6. Fix a map </?GOg(V, M) and assume that G(v o) is a non-degenerate critical 
G-orbit of p such that m~(C(v o)) ^ 0. Then for all open, G-invariant neighbourhoods 
0 C V of the orbit G(v 0 ) such that (V93) -1 (0)n@ = G(v 0 ) there exist an open, G-invariant 
neighbourhood U C cl(U) C 0 of the orbit G(v o), £ > 0 and a map <f> G Og(V, M) such 
that 

(1) p(v) = 4>{v) for all v G V\ U(e), where U(e) denotes an e-neighbourhood of the set 
U, that is U veu B v(V)> 

(2) G(v o) is a special critical G-orbit of d, 

(3) ((V(/>)- 1 (0) n (U(e)\G(vo))) <H) = (V^-^O) n (U(e)\G(v 0 )), where (. H) = (G Vo ). 
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Lemma 2.7. Fix a map <p€G/;(V, M) and assume that G(v o) is a non-degenerate critical 
G-orbit of tp such that m~(C(v o)) ^ 0 and (H) = (G Vo ). Then there exists an open, G- 
invariant neighbourhood 0 of the orbit G(v o) such that (V(/?) _1 (0) D0 = G(v o). Moreover, 

V G -deg(V<p,0) = (-l) m -W*>» X G(G/H + ) + 

+ ^ V G -deg (i? ,)(V</?, 0) • Xg(G/H' + ). 

(H')£sub{G],(H')<(H) 

Proof. We can choose an open, G-invariant neighbourhood 0 of the G-orbit G(uo) such 
that (Vy?) _1 (0)n@ = G(v 0 ). Without loss of generality we can assume (V<p) _1 (0)n<9@ = 0. 
Now, using Lemma 2.6, we define a ©-admissible homotopy H(v,t) = tp(v) + (1 — t)<f>{y). 
Next, the homotopy invariance property of the degree (Theorem 2.2.(3)) implies that 

V G -deg(Vy?, 0) = V G -deg(V„tf(-, 1), 0) = V G -deg(V t ,tf (•, 0), 0) = V G -deg(V0, 0). 

In addition, (V0)^ 1 (O) f!0 C U (|) and (V</>) -1 (0) fll/ (|) = {G(n 0 )}, just like in Lemma 
2.6. Therefore using excision and additivity properties of the degree (Theorems 2.2.(4) 
and 2.2.(2)) and Lemma 2.5 we can conclude that 

V G -deg(V0, 0) = V G -deg(V</>, U (e/2)) = V G -deg(V0, U (e/4))+ 

+ V G -deg(V0,C/(e/2)\cl(C/(e/4))) = (-1 )™P b ^)) Xg (G/ H + )+ 

+ V G -deg ( ^)(V0, U (e/2) \cl (U (e/4))) • X g(G/H’ + ), 

{H’)elub[G\,(H')<(H) 

where by cl (U (e/4)) we denote the closure of U (e/4). □ 

3. Statement of main results 

The following assumptions will be needed throughout this section. Let G be a compact 
Lie group and suppose that 12 C V is an open, G-invariant subset of a finite-dimensional, 
real, orthogonal G-representation V. Consider a G-invariant G 2 -potential p : 12 x M —y M. 
Since V v p is G-equivariant, if (v 0 , p 0 ) £ (V 7 ^^) —1 (0) then G(v 0 ) x {p 0 } C (V^<p) _1 (0). 
A set of the form G(v o) x {p 0 } consisting of critical points is called a critical G-orbit of 
ip. Using the techniques of equivariant bifurcation theory we will study solutions of the 
equation 

V v <p(v,p) = 0. (3.1) 

Additionally, assume there exists a continuous map w : M —> 12 such that 

P = [J G(w(p)) x {p} c (V^r'tO). 

peM 

The family T is called a trivial family of G-orbits of solutions of (3.1) and for any subset 
X C M let Tx denote the set of G-orbits [j peX G(w(p)) x {p} C T C 12 x 1. 

In this section we will formulate the necessary and sufficient conditions for the existence 
of bifurcations of nontrivial solutions of (3.1) from the trivial family T. 

Let us introduce a notion of a local bifurcation of G-orbits of solutions of (3.1). 
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Definition 3.1. Fix parameters p ± G M. such that p~ < p + . A local bifurcation from 
the segment of G-orbits T\ p - tP +] C T of solutions of (3.1) occurs if there exists a G-orbit 
T Po C J-[ P - )P +] such that the point (u;(po),po) G T Po is an accumulation point of the set 
{(u,p) G (12 x M)\J r : V„<p(u,p) = 0} (see Figure 1). 

We call po a parameter of local bifurcation and T Po a G-orbit of local bifurcation. We 
denote by EXT the set of all parameters of local bifurcation. 



Figure 1. A local bifurcation from the segment of G-orbits J-A- p +] C T. 

According to Lemma 2.1, for each parameter p G M. we have dimker V 2 <p(u;(p), p) > 
dim G(w(p)). Now we will formulate the necessary condition for the existence of parame¬ 
ters of local bifurcation. 

Theorem 3.1. Under the above assumptions, if, moreover, p 0 G EXT, then 

dimker V 2 v ip{w(p 0 ), p 0 ) > dim G(w(p 0 )). 

Proof. To prove this theorem we apply two kinds of implicit function theorems - the 
classical one and the G-invariant one. Suppose, contrary to our claim, that po £ EXT and 
dimker V,y(p(w(po),po) = dimG(tc(po)). Consider two cases: tc(po) G Y G and w(po) ^ V G . 
Case w(p Q ) G ¥ G . Since w(p 0 ) G V G , we have 

G(w(p 0 )) = {w(po)} and 0 = dim G(w(p 0 )) = dimker V^(p(w(p 0 ),Po), 

that is det V 2 <p(u;(po), Po) 7 ^ 0. The map <p is of class G 2 and w is continuous, so there 
exists e > 0 such that det V 2 (p(tc(p), p) 7 ^ 0 for all parameters p G (po — GPo + £)■ 
The implicit function theorem applied at point (w(pf),pf) yields that there exist open 
neighbourhoods U w ^ po ) C 12 and U Po C M. of the points u>(po) G 12 and po G M, respectively, 
and exactly one continuous map -0 : U Po —> U w ^ po ) such that 

V„<p(u,p) = 0 and (v,p) G G„, (po) x U po if and only if v = if(p). 


















GLOBAL BIFURCATION 


13 


Therefore w(p) G V G for all parameters p G U po , ip — w 


and hence 


0 


(V„<p) (o) n (U w{po) x U P0 \X) = 0, 


which contradicts that p 0 G BXT. 

Case w(po ) ^ V G . As a consequence of the G-invariant implicit function theorem (see [8]) 
and Remark 4 of [8] there exist an open, G- invariant neighbourhood 0 C of the G-orbit 
G(w(po)) and a continuous, G-equivariant map ip : G(w(p o)) x (p 0 — e, po + e) — > V such 
that if 

U = {(ip(v, p), p) : P G (p 0 - £, Po + e), U e G(w(p 0 ))}, 
then for all (v, p) G 0 x (p 0 — e, p 0 + e) 

V v <p(v,p ) = 0 if and only if (v,p) G A/”. (3.2) 


ffence for all p G (po — £, Po + s) if V.„<p(iq, p) = 0 and V„<p(u 2 , p) = 0, where tq, u 2 G 0, 
then tq = ip(v[,p) and v 2 = ip(v 2 ,p) for some i4 G G(w(po)). Thus there exists 
5 G G such that • Consequently, tq = ip(v[, p) = ip(gv 2 , p) = gip(v 2 , p) = gv 2 . 

On the other hand, we have the map w : M —>• Ql such that V v <p(w(p), p) = 0 for all 
parameters p G R. Thus for all p G (po - £,po + s) if w(p) G O then (■ w(p),p ) G A/”, 
that is w(p) = ip(gw(p 0 ), p) for some g G G. Since in is continuous and in(po) G 0, 
there exists small enough e < e such that w(p) G O for all p G (po — £, Po + £■)■ Hence 
(V„<p) _1 (0) fl (0 x (p 0 — i, po + e)\X) = 0, which again contradicts that po G £>ZX1 □ 


After introducing the local bifurcation of G-orbits of solutions of (3.1) we now turn 
to the notion of global bifurcation of these solutions. Denote by G(po) the connected 
component of the set cl({(v,p) G(ff x M)\X : V t ,<p(u,p) = 0}) U Xy J0 containing J r po and 
by G([p _ ,p + ]) the connected component of the set cl({(v,p) G (D x M)\J r : V„<p(v,p) = 
0}) UX[ p - iP +] containing X[ p -. p +]. 

Definition 3.2. Fix parameters p^ G M such that p~ < p + . A global bifurcation from 
the segment of G-orbits XA-^+i C X of solutions of (3.1) occurs if the component 
G([p _ ,p + ]) cSlxR satisfies the following condition: 

G([p _ ,p + ]) is not compact or (G([p“, P + ])\X[ p - iP +j) fl X^0 (3.3) 

(see Figure 2). 

We call a parameter p 0 G [p _ ,p + ] a parameter of global bifurcation if the component 
G(po) CfixMis not compact or (G(po)\Xj p - iP +]) fl X ^ 0. A G-orbit X po C Xj p - >p +] is 
called a G-orbit of global bifurcation. We denote by QCOB the set of all parameters of 
global bifurcation. 


The condition (3.3) is referred to as the Rabinowitz type alternative ([31]). Notice that 
QCOB C BXT. Below we define the bifurcation index from the segment [p , p + ] to be an 
element BIF^- p +] G U(G). 

Definition 3.3. Fix parameters p ± G M, p - < p + such that p ± ^ BXT and fix open, 
G-invariant neighbourhoods Q± C fl of G-orbits G(w(p ± )) such that 

(V„<p) _1 (0) n (cl(0 ± ) x {p ± }) = X p ±. 
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Figure 2. A global bifurcation from the segment of G-orbits JF\p- tP +] C T . 

We define the bifurcation index from the segment [p~,p + ] by 

BIF [p -, p+] = V G -deg(V t; <p(-,p + ),0 + ) - V G -deg(V^(-,p-),0-) G U(G). 

The definition of BIF\ p - jP +] does not depend on the choice of the G-invariant neighbour¬ 
hood of the G-orbit. In what follows we formulate a sufficient condition for the existence 
of global bifurcation from the segment of solutions of (3.1). 

Theorem 3.2. Fix parameters p ± G M, p~ < p + such that p ± (f BXF and assume 
that BIF[p-, p +] 7 ^ 0 G 17(G). Then a global bifurcation from the segment of G-orbits 
B \ p - iP +] C T of solutions of (3.1) occurs. 

Proof. First, we will prove that there exists a parameter p 0 G (p _ ,p + ) such that p 0 G BXF. 
Suppose, contrary to our claim, that BIF[ p - tP +j ^ 0 and BXF D [p _ ,p + ] = 0. Hence there 
exists an open, bounded G-invariant neighbourhood 7 C O x 1 of the set F\ p - jP +-\ such 
that (V l ,<p)^ 1 (0) D (c\(U)\F) = 0. In particular, we can pick open, bounded G-invariant 
neighbourhoods C H of G-orbits G(w(p ± )) such that (Vj,<p) - 1 (0) n(cl(@ ± ) x {p 1 * 1 }) = 
F p ±. The generalised homotopy invariance property of the degree (Theorem 2.2.(3)) yields 
that 

V G -deg(Vv?(-,p-),I/n (ft x {p-})) = V G -deg(Vp(-,p + ),Gn (fi x {p + })), 

while from excision property (Theorem 2.2.(4)) we obtain the following equalities: 

V G -deg(V<p(-,p-),t/n (fi x {p-})) = V G -deg(V<p(-,p-),0-), 

V G -deg(Vp(-, p + ), U n(Cl x {p+})) = V G -deg(Vp(-,p + ),0+), 

a contradiction. We have thus proved that there exists a G-orbit of local bifurcation F Po 
of solutions of (3.1), or equivalently, the parameter p 0 is a parameter of local bifurcation, 
that is p 0 G BXF. 

Secondly, we will prove that a global bifurcation from F [ p - jP +] occurs. The proof will be 
divided into two steps. 
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Step 1. We first prove that C([p , p + ]) ^ F \ P - )P +]■ Suppose the contrary and choose an 
open, bounded set Q C flxM such that J-) p - p +) C Q and fix ((—oo, p _ )U(p + , +oo))flQ = 
0. Additionally, let (fl x { y o =l= }) fl<9 Q = cl(0 ± ) x {p ± }, where 0 ± C fl are chosen similarly 
as above, that is 0“ C hi are open, bounded, G-invariant neighbourhoods of the G-orbits 
G(w(p ± )) such that (Vi,<p) _1 (0)n(cl(@ =b ) x {p ± }) = T p ±. Observe that cl(<5)n(V t ,<^)^ 1 (0) 
is compact and put 

K = cl(Q) 0 (V u <p) _1 (0), 

A = F\p -, p +], 

B = dQ D cl({(u, p) G (hi x M)\J r : V v (p(v , p) = 0}). 

These sets satisfy assumptions of Theorem 2.3, hence we obtain sets Ka and Kb from the 
conclusion of this theorem. Since these sets are compact and disjoint, there exists ■// > 0 
such that //-neighbourhoods Ka(p), Kb{t)) of the sets Ka, Kb are disjoint. Now put 

U = G(Q\cl(K B {^))) 

and observe that dU fl (V„<^) _1 (0) = T p ± by Lemma 2.4. The set U is G-invariant and 
open by Lemma 2.3. Again, using the generalised homotopy invariance property of the 
degree (Theorem 2.2.(3)), we obtain the equality 

V G -deg(V^(-,p-),0-) = V G -deg(Vv?(-,p + ),0 + ), 

a contradiction. 

Step 2. Now we prove that either the component C([p~,p + ]) is not compact or that 
the sets G([p _ ,p + ])\J r [ p - iP +] and T have nonempty intersection. Suppose the contrary 
and pick an open, bounded e-neighbourhood (cl(G([p _ , p + ])\J-'[ p - p +]))(e) C fl x M of 
cl(G([p-,p + ])\J r [p - iP+] ) such that cl((cl(G([p~, P + ])V>-,/>+]))( £ )) 0 ^(-oo, P -]u[ P +,+oo) = 0 
for some (small enough) e > 0. Let Q C fl x M and 0^ C fl be the sets chosen as in 
Step 1 satisfying an additional condition: (0“ x {p - })fl(cl(G([p _ ,p + ])\J-’[ p - iP +]))(e) = 0. 
Now define Q\ = Q U (cl(G([p _ , p + ])\J r [ p - iP +]))(e) and notice that cl(Qi) fl (V t ,p , ) _1 (0) is 
a compact set. Again put 

K = cltQi) n (V„<p) -1 (0), 

A = C([p~,p + ]), 

B = dQi nd({(w,p) G (fl x Mjyj 7 : V„<p(u,p) = 0}). 

The rest follows analogously to Step 1, but with Q\ instead of Q. So defined sets satisfy 
the assumptions of Theorem 2.3, therefore we obtain Ka and Kb in an identical manner as 
earlier. Once again they are compact and disjoint, which enables us to pick their disjoint 
//-neighbourhoods Ka{t)) and Kb{t])- Similarly to Step 1, we now put 

C/ = G(Q 1 \d(A' B (t I ,))) 

and notice that U is open, G-invariant and dU fl (V„<p) -1 (0) = T p ±. Once again the 
generalised homotopy invariance property of the degree (Theorem 2.2. (3)) implies equality 
of degrees, a contradiction. □ 
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A consequence of Theorem 3.1 is that there is no G-orbit of any type of bifurcation 
among non-degenerate G-orbits, which is made precise by the following: 

Corollary 3.1. Fix parameters p ± ,Po £ M, p _ < po < p + and assume that T f>0 is the 
only degenerate G-orbit in the segment [p~, p + ]. Thus, if BI F[ p - )/0 +] / 0 6 17(G), then 
there exists exactly one G-orbit of local bifurcation in the segment [p~,p + ] and this orbit 
is precisely J : PQ . Moreover, it is also a G-orbit of global bifurcation. 

What remains is a question how to distinguish the degrees of critical orbits. It transpires 
that under a very mild assumption this can be done by means of a simple condition. 

Theorem 3.3. Fix parameters p ± £ M, p~ < p + and assume additionally that the G- 
orbits G(w(p ± )) are non-degenerate and (G w ( p - iq) = (G w ( p -)) = (77). Thus, if 
det B{w(p~)) ■ det B(w(p + )) <0, then 

V G -deg(V„<p(-,p-),©-) ± V G -deg(Vpp(-,p + ),0 + ). 

Proof. According to Lemma 2.1, there is a special form of V^(p(in(p =b ),p^) given by the 
formula (2.3). Observe that if m~ (G(w(p ± ))) = 0, then 

V G -deg(V^(-,p ± ),0 ± ) = (:_i)--(b(-(p ± ))) Xg( g/77+). (3.4) 

When, on the other hand, m _ (G(w(p ± ))) ^ 0, by Lemma 2.7, we obtain that 

v G -deg (v^t-.pq.eq = (-ir-m^ ± m Xa{G/H + ) + 

+ 5] Vc-deg^jtV^f., p), e±) • Xg(G/H'+). 

(3.5) 

Since det B(w (p - )) ■ det B(w (p + )) < 0, we get that (—l) m ( b G’(p ))) ^ (—l) m ( b ( w (p + )))_ 
Using the formulas (3.4) and (3.5) we complete the proof. □ 

Remark 3.1. In the case of Theorem 3.3 with (G w ( p -)) ^ (■ G w ( p +)), we get inequality of 
the degrees V G -deg(V t ,<p(-, p ± ), 0 1 * 1 ), what is the consequence of the formulas (3.4) and 
(3.5). 

The following theorem is in the spirit of Rabinowitz’s theorems regarding the problem 
of global bifurcation of zeros without symmetries (see [31]). 

Theorem 3.4. Fix parameters p ± €M, p~ < p + such that p ± BTT and BIF\ p - p +] ^ 
0 £ 17(G). Then a global bifurcation from the segment of G-orbits 7 7 \ P -, P +] C P of so¬ 
lutions of (3.1) occurs, that is the component G([p“,p + ]) C 12 x M is not com.pact or 
(G([p _ , P + ])\7 r [p- )P +]) n J / 0. Moreover, if the component G([p“,p + ]) is compact and 
there exist finitely many parameters 

Po < Po < • • • < Pi = P~ < P + = pt <■■■< P~ < Pn 

such that 

(1) C([p-, G ])n^cULop,4.^i. 

( 2 ) C{\p-,p])nr Kjt \ # 0 for k = 0 , 
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(3) GCOB n ULo-F^+i = cl(C([p-,p + ])\^j,,-^]) n ULopp; ,/>/]• 

(4) BIT for k = 0,..., n, 
then 

n 

0 - 6 ) 

fc=0 

Proof. According to Theorem 3.2, we only need to show that the formula (3.6) is valid. 
To prove this, fix an open and bounded ^-neighbourhood (cl(G([p _ , p + ])\-F[p-,p+] ))(£), de¬ 
noted by W for short, of the set cl(G([p _ , p + ])\-F[ /3 -,p+]) such that cl(4U)nJ-’ R y|jri ^ p - p +^ = 
0 for some (small enough) e > 0. By assumption (3), for any k 4 l we can choose an open, 
bounded £-neighbour hood (cl(G([p)7, p +] UC([p“,p + ]))))(e), denoted by 114 for 

short, of cl(C([pfc , Pfc ])\(-? r [p-, p + ] UC([p", p + ]))) such that cl(W fc ) n^R\ ( p-, p +) = 0 for some 

(small enough) e > 0. For any k G {0,..., n} set an open and bounded set Qt C O x R 

such that J 7 ( p -. p +) C Qk and (hi x ((— oo,p k ) U (p£,+oo))) D Qk = 0- Additionally, let 

(O x {p^}) fl dQ k = cl(0 p ) x {p p }, where 0 P C hi are an open, bounded, G-invariant 

+ .. 

neighbourhoods of the G-orbits G(w(p k )) such that (V u <p) -1 (0) fl (cl(0jr) x {p k }) = 
T p ± and (0^ x {p^}) fl W = 0. For k ^ l let 0 P satisfy an additional condition 

(®fc x 0 W'fe = 0- Now for k ^ l dehne Qk = Qfc U 144 and for k — l let Qi = Qi. 

Next, dehne Q = [j k=0 Qk U W and observe that cl(Q) fl (V t ,<p)^ 1 (0) is a compact set. 
Put 

A' = ci(g)n(v u <p)~ 1 (o), 

n 

^=U^k,* uC ([X,p + ]) 

k =0 

B = dQ fl cl({(v, p) G (U x M)\A : V„<p(u, p) = 0}). 

Sets so dehned satisfy the assumptions of Theorem 2.3, hence we obtain sets Ka and 
Kb from the conclusion of this theorem. Since these sets are compact and disjoint, there 
exists rj > 0 such that //-neighbourhoods K^(rj), Kb{t]) of the sets Ka , K B are disjoint. 
Now put U = G(Q\c\(K B (lrj))) and observe that (V„<p) _1 (0) fl dU = IJ!L 0 J 7 - 0 ' + 
by Lemma 2.4. The set U is G-invariant and open by Lemma 2.3. Using the generalised 
homotopy invariance property of the degree (Theorem 2.2.(3)) we obtain that 

V G -deg(Vpp(-,p),Un (fi x {p})) G U(G) 

is well dehned and is constant as a function of the p G [p Tmn + e 0) Pmax ~ £o]> where 

Pmin = inf{p Gl:h n(Ox {p}) ^ 0}, 

Pmax = SUp{p GlT 0(OX {p}) ± 0} 

and £o is a positive number satisfying (V„<p) -1 (0) fl (C/fl (U x {p m in ±£o})) = 0 hi the case 

max 

of Pmin < Po and Pmax > p£. In Other Cases, p G [P(7, Pmax ~ £q] ° r P £ [pmin + £q, P%] Or 
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p G \p 0 , Pn]- Observe that, using the additivity property of the degree (Theorem 2.2.(2)), 
we have 

V G -deg(V^(-,p+),f/n(Ox{p+})) = 

V G -deg(V„<p(-,p+), 0+) + V G -deg(Vpp(-,p+), (17 n(Ox {p+}))\cl(0+)), 

V --v-' - V -- — — * 

a t C k 

analogously for parameters p^, k — 0,..., n. We assume that V G -deg(V„<p(-, p), 0) = 0 G 
U{G). Since = a)7 +1 + c^ +1 , we have a~£ = % +1 , hence we conclude that 

n n n n n— 1 n— 1 

a k ~ a k = a k — = °fc+l = — a fc+l ) + a n ~ a 0 = °n — °0 ’ 

fc=0 fc=0 fc=0 k =0 fc=—1 fc=0 

Observe also that c+ = Cq = 0. Then 


£ = £ v G .-deg(v,,y>(., P +),e+) - v G -deg(v„p(., Pl -),e l :) = 

k =0 k =0 

n 

= _ = a n — a 0 = ( a n + C n) ~ (°0 + C 0 ) = 0 G U ('G), 

fc =0 

which proves our assertion. □ 

Remark 3.2. Observe that in Theorem 3.4 if for k — 0,..., n the set 7L- „+i contains 
exactly one degenerate critical orbit, then the assumption (3) of this theorem is satisfied. 

Now fix (v 0 , Po) G (V.„<p) _1 (0) and let 

V+ = (J) V_vX„ 0 ,p 0 )(A), 

AS(T+(-V 2 ¥ 3 (?Jo,Po)) 

where V-v2<p(v 0 ,po)W denotes the eigenspace of — V 2 tp(u 0 ,Po) associated to A. Notice that 
a non-degenerate critical G-orbit G(v o) C fl of the map <p(-, p 0 ) is an isolated r/ po -invariant 
set, where r/ Po : U C fi x R G ¥ is a local G-flow induced by the equation 

v(t) = —V v <p(v(t),p 0 ). 

For simplicity of notation, G Vo = H. According to Lemma 2.1, there is a special form of 
V 2 <p(uo,po) given by the formula (2.3), where 

V = T V0 G(v 0 ) © W H © (W^) 1 - = T V0 G(v 0 ) © V+ © (V + ) ± . 

Assume that the non-degenerate critical G-orbit G(v 0 ) is special, that is m~(C(v 0 )) = 0. 
Then the G-equivariant Conley index of the critical G-orbit G(v o) has the following G- 
homotopy type: 

«c(GW,1„) = ([(G x f; D(V+))/(G x„S(V + ))] G ,[*]), 
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where D(V + ), S(¥ + ) denote a closed ball and a sphere in the space V + , respectively. Ob¬ 
serve that since m~(C(v 0 )) = 0, we have S(Y + ) C D(Y + ) C V + C W H and m~(B(v 0 )) = 
dimV + , hence 

CWm),?,) = ([{G/H x D(V+))/(G/H x 5(V + ))]o, [»]). (3.7) 

Remark 3.3. Additionally, there is a connection between the G'-equivariant Conley index 
of the special critical G-orbit G(v o) and a notion of finite, pointed G-CW-co rnplex. The 
G-equivariant Conley index CIg(G(vo),t] P o ) has the G-homotopy type of a pointed G- 
C IT-complex which consists of the base point * and one rn~ (R(n 0 ) )-dimensional cell of 
orbit type (H). 

Remark 3.4. Thus, if Morse indices of two special critical G-orbits are different, then 
the G-equivariant Conley indices have different G-homotopy types. It is known (see [37]) 
that a change in the G-equivariant Conley index implies the existence of G-orbits of local 
bifurcation of critical points of G-equivariant variational problems. 

As a consequence of Remark 3.3 and Remark 3.4, there are simple conditions verifying 
the existence of a local bifurcation. 

Theorem 3.5. Fix parameters p ± GM, p~ < p + such that: 

(1) dimker 'Vl(p(w(p ± ), p±) = dim G(w(p ± )), 

(2) (G w ( p -)) = (G W ( P +)), 

(3) m~(C(w(p ± ))) = 0. 

If m~ (B(w(p~))) ^ m~ (B(w(p + ))), then (p ~, p + ) FEXB ^ 0, that is there exists a local 
bifurcation parameter p 0 G (p~, p + ). 

Proof. First observe that, according to Remark 3.3, the G-equivariant Conley indices of 
the special critical G-orbits G(tc(p =b )) have the G-homotopy type of a pointed G-CW- 
complex which consists of the base point * and one rrr ( B (w(p ± )) )-dimensional cell of 
orbit type (G w ( p -)) = (G w ( p +)). Since m~(B(w(p~))) ^ m~(B(w(p + ))), we conclude, by 
Remark 3.4, that (p~, p + ) D BXIF ^ 0. □ 

Theorem 3.6. Fix parameters p ± GM, p~ < p + such that: 

(1) dimker V 2 v (p(w (p ± ),p ± ) = dimG(w(p ± )), 

(2) (G w ( p -)) 7^ (G M (p+j), 

(3) m~(C(w(p ± ))) = 0. 

Then (p~, p + ) D EXT ^ 0, that is there exists a local bifurcation parameter po G (p~, p + ). 

Proof. Analogously to the proof of Theorem 3.5, the G-equivariant Conley indices of the 
special critical G-orbits G(w(p ± )) have the G-homotopy type of a pointed G-GIT-complex 
which consists of the base point * and one m“(R(w(p ± )))-dimensional cell of orbit type 
(G w ( p ±)), see Remark 3.3. Since (G w ( p -)) ^ (G w ( p +f), applying Remark 3.4 once again, 
we get that {p ~, p + ) D BXB ^ 0. □ 
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4. Applications 


In this section we apply our results to show the existence of new families of central 
configurations, which bifurcate from certain known families of central configurations. We 
consider N bodies of positive masses mi,..., mjv in the 2-dimensional Euclidean vector 
space, whose positions are denoted by q±, G M 2 . The equations of motion are given 
by 


m/lj 


N 

^ G mjrrii 


Qj ~ Qi 


j = 1,..., N, 


where G is the gravitational constant, which can be set equal to one. From now on we 
will treat the space M. 2N as an 50 (2)-representation V, that is a representation of the 
Lie group 50(2), which is a direct sum of N copies of the natural, orthogonal 50(2)- 
representation (that is the usual multiplication of a vector by a matrix). The action of 
50(2) on V is given by 


50(2) xY 3 (g,q) = (g, (gi,..., q N )) ^ g ■ q = (g -ft,..., g ■ q N ) e Y. 


Newtonian equations of motion associated with a potential U have the following form: 


Mq = W q U(q,m), 


(4.1) 


where M is a mass matrix, that is M = diag(rrii,mi ,...,mjv,^v), and the Newtonian 
potential U : fl x (0, + 00 )^ —> M is given by 


U(q,m) = U(q 1 ,.. .,q N , m 1 ,.. ., m N ) 


£ 

1 <i<j<N 


mirrij 

I Qi -qjY 


where the set C V is defined by 


^ = {?= (9i, • • • \ Qn) e ¥ : qi ± q 3 for i ± j}. 

We can assume without loss of generality that the centre of mass of the bodies is at the 
origin of the coordinate chart. Observe that the set 0 C V is 50(2)-invariant and open, 
while the potential U is an 50(2)-invariant 0°°-map. 


Definition 4.1. A configuration q — (qi,, q n) G 0 is called a central configuration of 
the system (4.1) if there exists a positive constant A such that q = —Xq. 


This is equivalent to saying that the following condition is fulfilled: 

- W q I(q,m) — V q U(q,m), (4.2) 

where I : fl x ( 0 , +oo) N —> M given by the formula 

1 N 

I{Ym) = 2 X!' ( ‘i 2 

3 =1 

is the moment of inertia. Using Euler’s theorem and the fact that U is homogeneous of 
degree —1 one can show that A = Moreover, for any central configuration q it is 

known that rq and g-q are central configurations with new A coefficients equal to ^3 and A 
respectively, where ret and g G 50(2). Thus this way we can introduce an equivalence 
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relation by saying that two central configurations are equivalent if we can pass from one 
to another by a composition of a rotation and a scaling. However, instead of passing to 
the quotient and treating classes as single points we shall study the whole 5'0(2)-orbits of 
central configurations. In this approach our problem becomes a problem of finding critical 
S'0(2)-orbits of a smooth, SO (2)- invariant potential <p : fl x (0, +oo) N —> M given by 

<p(g, m) = U ( q , m) + A J(g, m ). 

As we mentioned earlier we search bifurcations of central configurations from known 
families, so we assume additionally that there exist continuous maps w : (0, +oo) —>■ fl 
and m : (0, +oo) —> (0, + 00 )^ such that 

V q <p(w(p),m(p)) = 0 . 

Now define a potential cp : D x (0, + 00 ) —>■ M by 

<p(?,p) = <p(9,™(p))- 


Then 

T= U SO(l)(«W)xWcW 1 (0). 

pe(0,+oo) 

We shall refer to T as the trivial family of solutions. Thus, we will consider the following 
equation: 

V g ^(g,p)= 0, (4.3) 

where A = A(p) = ■ Now we will apply our main Theorems 3.1 and 3.2 to 

formulate necessary and sufficient conditions for the existence of bifurcations of central 
configurations from our known families w. Note that for families w considered in this paper 
we have SO( 2 )„,( p ) = Id thus the matrix C[w{p )) is zero dimensional and in particular, 
its Morse index is equal to zero. Now, as has been indicated earlier, we pass from seeking 
bifurcations of central configurations to seeking bifurcations of critical S'0(2)-orbits of the 
»S'0 (2)-invariant potential 

< P( < ?) P) = U(q, m{p)) + A (p)I(q,m(p)). 


For this we seek parameters p satisfying the necessary condition for the existence of a 
local bifurcation, which is given by Theorem 3.1. Because of tedious computations we 
aid ourselves by the algebraic processor MAPLE 1M to calculate the Hessian V ^<p of the 
potential <p and its characteristic polynomial W p along the curve w. Denote the symbolic 
form of the latter by 

W p (q) = q 2N - a\(p)q 2N ~ l H-f (- l) 2N ~ 1 a 2N -i(p)q + (~1) 2N a 2N (p)■ 

Notice that dimker V 2 ip(w(p), p) = k if and only if a 2 N = ••• = a 2 N-k+i = 0 and 
a 2 N~k 7 ^ 0 . Additionally in this case a 2 M~k is the product of the nonzero eigenvalues. Also 
recall that dimker V 2 ip(w (p),p) > dim SO(2)(w(p)) = 1, so a 2 7 v(p) = 0. For this reason 
we consider the equation 

fl2A-i(p) = 0, 

whose solutions are the only parameters satisfying the necessary condition from Theorem 
3.1. In the next step we apply Theorems 3.2, 3.3 and 3.5 to each candidate for the 
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parameter of bifurcation to separate the ones which also satisfy the sufficient conditions. 
For each candidate p* we choose its small enough e-neighbourhood for which we have 
dimker V 2 q q)(w(p* ± e),p* ± e) = dim SO(2)(w(p* ± e)). If 

m~(B(w(p* - e))) ^ m~(B(w(p* + e))), 

then there exists a parameter of local bifurcation in this neighbourhood. Moreover, if 
these Morse indices differ by an odd number a global bifurcation occurs. 



Figure 3. Eight bodies located at the vertices of two nested squares. 


We will now focus on our goal, which is studying bifurcations from the known families 
of central configurations, which we will now consider as our trivial family of solutions w. 
First, we consider nested planar central configuration for the problem of 8 bodies. This 
configuration was found in 2013 by A. C. Fernandes, L. F. Mello and M. M. da Silva 
in [12] and consists of 4 bodies with equal masses located at the vertices of a square 
whose side is equal to 1 and the other 4 bodies, also with equal masses, located at the 
vertices of a smaller square whose side is equal to 0 < b < 0.53177..., whose centre 
coincides with the centre of the first square (see Figure 3). Denote by r the ratio of the 
circumcircle of the inner square. We treat 0 < r < tq = 0.37602... as a parameter. Let 
mi = m 2 = m 3 = m .4 = M r = — W(rq anc ^ 1715 = m 6 = 1717 = 757,8 = ^ where 


A(r) 

B(r) 


(R\,2 — 7?i ! 5)Ai i 5 i 2 + (Ri,3 ~ -Rl, 6 )^ 1 , 6 , 3 + (Rl,7 — -Rig)^1,7,2, 
(Re, 7 — Rl, 5 )^ 1 ,5,6 + (Rl,7 — Re, 7 )^ 5 ,6,3 + (Rl,6 — R$, 7 )^ 1 ,6,8 


and R id = 1/|^ - qj | 3 , A iJtk = (g, - qj) A (qt - q k ) for 1 < i,j,k < 8. We define 
w : (0, r 0 ) -A D by 


...111 111 11 y/2 V2 V2 Vz V2\ 

w[r)=\ —, —- r, -r, r,-r,-r, r,-r, r 

' ' i 2 ’ 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2/ 

(4.4) 
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and m : ( 0 , r 0 ) —» ( 0 , +oo ) 8 as follows 

m(r) = (M r , M r , M r , M r , 1,1,1,1). 

Lemma 4.1. Put r\ = l\/2j u 2 = |\/2 and r 3 = |\/2. T/ien dimker V 2 </?(u;(rj), rf) = 
dim 5'0(2)(in(rj)) = 1 for i = 1,2,3 and the Morse index of V 2 tp(w(-), •) evaluated at r t 
is 

( 1, for i — 1 

m_ (V^(w(rj),ri)) = < 3, for i = 2 . 

[ 4, for f = 3 

So, there is no matrix C(w(r)) and 50(2) w ( r ) = Id for any 0 < r < r 0 = 0.37602... 
Also, by Lemma 4.1, we have dimker V 2 </?(in(?y),rj) = dimS'0(2)(ta(rj)) = 1 for i = 
1,2,3. Thus, by Theorem 3.5, since m~(V q <p(w(ri), rr)) 7 ^ m~(f\7 2 q (p{w{rf),r2))-, there 
exists a local bifurcation parameter in the segment (r 1 , r 2 ). Additionally, by Theorem 3.2 
and Theorem 3.3, since the numbers m _ (V 2 </?(u;(r 2 ), r 2 )) and m _ (V 2 </?(w(r 3 ), r 3 )) are of 
different parity, we have that there exists a global bifurcation parameter in (r 2 , r' 3 ). 

Summarising, we have proved the following theorem: 

Theorem 4.1. 

(1) There exists a local bifurcation parameter in the segment (rr,r 2 ), that is (rr,r 2 ) D 

7^ 0. 

(2) There exists a global bifurcation parameter in the segment (r 2 ,r 3 ), that is (r 2 , r 3 ) D 
^£0/5 7^ 0. 

Notice that we can consider a subset of the full configuration space 12 which is in¬ 
variant for the gradient flow (the set of central configurations of two nested squares type 
(r, m, M)), then studying central configurations in this set becomes a problem of studying 
zeros of a function F : (0, r 0 ) x (0, + 00) 2 —> M given by the formula F(r, m, M) = MA(r) + 
mB{r). For the trivial family of solutions (r, m, ), for any (r, m) G (0, r 0 ) x (0, + 00 ), 

we have F (r,m, M(r,m)) = 0, where M(r,m ) = and F' (r, m, M (r, m)) = 

M(r,m)A'(r) + mB'(r) > 0, so by the implicit function theorem there is no bifurcation 
of central configurations of two nested squares type from the trivial family. Combining 
Theorem 4.1 and the above, we obtain that the families which bifurcate from the family 
(4.4) are not of two nested squares type. 

Now consider a rosette central configuration consisting of n particles of mass m 1 lying 
at the vertices of a regular n-gon, n particles of mass m 2 lying at the vertices of another 
n-gon, where the second one is rotated by -, and an additional particle of mass mo lying 
at the common centre of the two n-gons (see [19] and [34]). We analyse bifurcations from 
this family in the case of 13 bodies (see Figure 4) which are considered in the following 
positions: 

q k+ i = <F (jfk) (n, 0), q 7+k = <E> (§ + ?fk) (r 2 , 0) for k = 0,..., 5 and q i3 = (0, 0), 
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7 ) 1 , 


7 ) 1 , 


1 ( 



11 


7 ) 1 . 


n i. 


12 


Figure 4. Thirteen bodies in rosette configuration. 


where by we denote the universal covering of SO (2) by M, that is a mapping $ : M —> 
SO (2) given by the formula 



We transform the coordinates rr and r 2 in the following way: rr = r cos 9, r 2 = rsin#, 
where ( r,9 ) G (0, +oo) x (0, |) (see [19]). We choose 6 = |, r = 1 and define w : 
( 0 , +oo ) 2 —> hi by 


m(m 0 ,mi) = (?i,..., <j 12 , q I3 ), 


(4.5) 


where masses mo and m\ are treated as parameters. Notice that for the two-parameter 
family w we obtain the same results as in Theorem 3.1 and Theorem 3.2. Namely, for any 
points and in the parameter space we can choose a one-parameter 

path joining these points, that is the segment connecting them. Then the configuration 
q— ((ft,..., qu, q 13 ) is central for each parameter (m 0 , mi) G (0, +00) 2 if 


1 


m 2 = 7n 2 (m 0 ,m 1 ) 



1862V3 - 7203 + 810^ + 90^^ 


So, we define m : (0, + 00) 2 —> (0, +cxd) 13 by the formula 

m(m 0 , mi) = (mi, mi, mi, mi, mi, mi, m 2 (m 0 , mi),m 2 (m 0 , mi),m 2 (m 0 , mi), 

m 2 (m 0 , mi), m 2 (m 0 , mi), m 2 (m 0 , mi), m 0 ) 

and compute the characteristic polynomial W( m0imi ) along the curve w. Then, we get that 

a 26 (m 0 ,mi) = 0 


and 


%(m 0 ,mi) = C(m 0 ,mi)C , 1 2 (mo,mi)C' 2 (mo,mi)C , |(m 0 ,mi), 
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where polynomials C'i, C 2 and C 3 are equal to 0 on some curves (see Figure 5) and 
C(mo,mi) 7 ^ 0 for G (0,+oo) 2 . 



- C l( m 0’ m l) =O 

« j )=0 

-Uj( m O’” , l) =0 


Figure 5. The set of zeros of the coefficient 025 and regions C3, C5, C2 and Co where 025 ^ 0. 


Lemma 4.2. The Morse index ofVfa along w depends on (mo, mj) as follows 


m (Vq<p(w(mo, mi), (m 0 , mi))) 


3, for (m 0 ,mi)ec 3 
5, for (m 0 ,mi)ec 5 
2 , for (m 0 ,mi)ec 2 
0 , /or (m 0 ,mi)ec 0 


and 


dimker V 2 <y9(ta(m 0 ,mi), (m 0 ,mi)) = dim5 , O(2)(ta(m 0 ,mj)) 
for (m 0 , mi) G (0, +oo) 2 such that Ci(m 0 , mi) 7 ^ 0 for all i = 1,2, 3. 


1 


So, for any (m 0 ,mi) G (0, +00) 2 we have that 5'0(2) U) ( mo!mi ) = Id and there is no 
matrix F(w(mo,mi)). Lemma 4.2 now shows that dimker V 2 9?(in(mo, mi), (mo, mi)) = 
dim < S'0(2)(iu(mo,mi)) = 1 for (m 0 , mi) G (0, +00) 2 such that Cj(m 0 , mi) 7^ 0 for all 
i G {1, 2, 3}. Thus, by Theorem 3.5, since the numbers m _ (V 2 99(in(m 0 , mi), (m 0 , mi))) are 
different in the regions c 3 , C5, C2 and Co, we have (mo, mf) G /3XJ 7 for (m 0 , mf) G (0, +00) 2 
with Cj(mo, m i) — 0, where i G {1,2,3}. Moreover, by Theorem 3.2 and Theorem 3.3, 
since the numbers m _ (V 2 </?(u;(mo, m 3 ), (m 0 ,mi))) are of different parity in the regions c 3 
and C2, we have (m 0 ,mi) G QCOB for (m 0 , mi) G (0,+oo ) 2 such that C2(mo,mi) = 0. 

Summarising, we have proved the following theorem: 


Theorem 4.2. 

(1) For any i G {1,2, 3}, if Ci(mo, mi) = 0, then (mo, mi) G BIJ 7 . 

(2) //C2(mo,mi) = 0, f/ien (mo, mi) G QCOB. 
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Notice that we can consider a subset of the full configuration space hi which is invariant 
for the gradient flow (the set of central configurations of rosette type mo, mi, m2)), 
then studying central configurations in this set becomes a problem of studying zeros 
of a function F : (0, +00) 3 — y M, F(x,£,fx) given by the formula (6) from [19], where 
x — ^ and n = For the trivial family of solutions (4.5), for any mo, mi G 

(0, +00), we have F (a/ 3, ma( ” ) , = 0 and F' (^3, > 0, so by the 

implicit function theorem there is no bifurcation of central configurations of rosette type 
from this family. Combining Theorem 4.2 and the above, we obtain that the families 
which bifurcate from the trivial one are not of rosette type. 



Figure 6. One-parameter family of rosette configuration where mi is a parameter and the 
set of four parameters of local bifurcation. 


Now we £x mi) G (0, + 00 ) and treat w(m^, •) and m(mj, •) as one-parameter families 
where rri\ is a parameter. Observe that there are exactly four parameters mi for which 
critical orbits SO(2)(w(mQ,mi)) are degenerate and there is at most one parameter m\ 
such that (mj,m[) G QCOB (see Figure 6). Then choose £ > 0 (small enough) such that 
(m^, m* ± e) ^ BXF and there is only one degenerate critical orbit SO(2)(w(mQ,m : [)) in 
the segment [m* — £, m\ + e]. It follows that 

BIF [ m *_ £)Tn * +e ] = VG-deg(V„<£>(•, (mo,m];+£)),0 + )-V G -deg(V I ,^(-, K,m[-e)), 0") = 
= ( _ 1)m -(^K^+-))) X50(2) ( > so(2)/{/d}+)-(-i)--(B(-K, m i- £ ))) Xso(2) ( >SO (2)/{/ rf }+) 

= ((-l) m “( S ^( m 5’ m l+e») _ (_i)m-(B(«;(m5,mI- C )))j Xgo(2 ^ jS - 0 ( 2 j^ /d |+j = 

= ((-l) 3 - {-lf)xso ( 2){SO{2)/{Id} + ) = -2xso(2)(SO(2)/{Id} + ) ? 0 G 17(G). 

Similarly, for any rn\ 7^ m\ since the numbers m _ (5(w(mj,mi ± d))) are of the same 
parity, we conclude that BIF[ mi _^ mi+ ^ = 0, where e is chosen as above. Therefore, by 
Theorem 3.4, C([m* — £,m^ + e]) is not compact. 
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